RichardsFoam3: A new version of RichardsFoam for continental surfaces hydrogeology modelling

RichardsFoam3 is an updated version of the OpenFOAM® solver RichardsFoam, previously presented in “An open source massively parallel for Richards equation: Mechanistic modelling water fluxes at watershed scale” by L. Orgogozo, N. Renon, C. Soulaine, F. Hénon, S.K. Tomer, D. Labat, O.S. Pokrovsky, M. Sekhar, R. Ababou, Quintard (Comput. Phys. Commun. 185 (2014) 3358-3371, https://doi.org/10.1016/j.cpc.2014.08.004), and new announcement “RichardsFOAM2: a RichardsFOAM devoted to vadose zone” Orgogozo 196 (2015) 619-620, https://doi.org/10.1016/j.cpc.2015.07.009). This includes improvements memory handling on-the-fly control computations, better integration framework, simplifications coding some expressions, as well advanced boundary conditions. All together these developments allow enhance ease application code continental surfaces hydrogeology modelling, its computational performances readability. The description elements contained this release may be found readMe file. Please note that you also find on hydrology page develop.openfoam.com interface: https://develop.openfoam.com/Community/hydrology/ Program title: CPC Library link program files: https://doi.org/10.17632/vkr7sd6fhb.1 Developer's repository link: Licensing provisions: GPLv3 Programming language: C++ Supplementary material: full system equations solved documentation directory package, with numbering used comments main file RichardsFoam3.C RichardsFoam3, order identify which lines are related equation. Journal reference previous version: Comput. 619-620 Does supersede version?: Yes Reasons Introducing functionalities implementing ergonomy, readability Summary revisions: features compared [1] RichardsFoam2 [2] following: (i) use during computation; (ii) possibility modification parameters linearisation procedure; (iii) conditions compiled along itself, dedicated simulation rain infiltration (rainFallFlux) or no (noFlux) while allowing exfiltration stand-alone (i.e. without mandatory swak4foam it was case [2]); (iv) Computation field piezometric head computation. NB: demonstration cases, way implement constant OpenFOAM proposed. (v) functions build seasonally variable easily; (vi) postprocessing procedures (vii) several other minor rewritings/cleanings code. Nature problem: software solves non-linear three-dimensional transient equation, very popular model transfer variably saturated porous media (e.g.: soils). It designed take advantage computing performance OpenFOAM®. goal able natural hydrosystems large temporal spatial scales. Solution method: A mixed implicit (FVM object oriented framework) explicit (FVC discretization equation backward time scheme coupled linearization method (Picard algorithm). Due loop final solution each step tends towards fully solution. implementation has been carried out concern robustness efficiency. Quintard, An mechanistic scale, 3358–3371, https://doi.org/10.1016/j.cpc.2014.08.004. RichardsFOAM2: zone, 619–620, https://doi.org/10.1016/j.cpc.2015.07.009.